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Abstract 


We consider the evolution of electromagnetic fields coupled to conduction currents during the 
reheating era after inflation, and prior to the establishing of the proton-electron plasma. We assume 
that the currents may be described by second order causal hydrodynamics. The resulting theory is 
not conformally invariant. The expansion of the Universe produces temperature gradients which 
couple to the current and generally oppose Ohmic dissipation. Although the effect is not strong, 
it suggests that the unfolding of hydrodynamic instabilities in these models may follow a different 
pattern than in first order theories, and even than in second order theories on non expanding 
backgrounds. 
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I. INTRODUCTION 


in 


The existence of magnetic fields on galactic and larger scales is one of the main puzzles 
present day cosmology [l-3|. Neither of the two major paradigms proposed to attack this 


question, namely dynamo amplification and primordial origin, seems to be able to provide 
a solution by itself jd]. It therefore seems likely that both mechanisms are at work, i.e., a 
seed field is generated early in the cosmic evolution and then subjected to one or several 
amplihcation st^es j^. This calls for a careful analysis of the cosmological history of 
magnetic fields |6|. 

Lots of efforts have been made to understand the evolution of primordial helds in the 
proton-electron plasma during the radiation dominated epoch. Special mention deserves the 
studies that address turbulent evolution, where helds with non-trivial topology i.e., with 
non zero magnetic helicitw would not be washed out by expansion as quickly as those with 
null magnetic helicity jb-S]. 

If we accept the existence of Inhation, then there must be a stage between it and the 
establishing of the proton-electron plasma where non-equilibrium processes dominated. That 
epoch is known as ‘reheating’. Moreover, electroweak (EW) and quantum-chromodynamic 
(QCD) phase transitions could have taken place by the end of it. Little is known of this 
epoch, besides the fact that all matter is created by the oscillatory decay of the inhaton. 
For example, the typical relaxation times and correlation times of the different interactions 
are not known. 

In this paper we shall perform a preliminary (see below) analysis of the evolution of 


magnetic helds during the reheating era 


10| . To this end, we shall consider that, on top 


of the two dominant contributions to the energy density, namely the coherent oscillations 
of the inhaton [u] and the incoherent radiation held, there is a charged huid that may 
interact non-trivially with the electromagnetic held. We do not identify this huid with the 
usual proton-electron plasma because we consider the evolution during an epoch well before 
quantum-chromodynamic phase transition. 

Both the coherent electromagnetic helds and the charged huid could be created as a 


side ehect of reheating by the parametric amp 
scalar held, as it has been discussed elsewhere 
could be the lightest supersymmetric partner, the s-t [12[. We shall assume that this huid 


ihcation of vacuum huctnations of a massive 
7|. A suitable candidate for the massive held 
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supports both viscous stresses and conduction currents, namely, electric currents without 
mass transport. For simplicity, we shall use Maxwell theory to describe the helds, in spite 
of the fact that the temperatures involved may be above the electroweak transition. 

At those early epochs the temperature and curvature of the Universe are very high and 
consequently a generally relativistic treatment is mandatory. The theory of relativistic real 
fluids has a long history but only relatively recently it has been put to the test, through its 
application to relativistic heavy ion collisions (RHICs) 13|]. Simply put, the straightforward 
covariant generalization of the Navier-Stokes equations leads to the so-called hrst order 
theories (FOTs), of which the Eckart 1^ and Landau-Lifshitz 1^ formulations are the best 
known. These theories have severe formal problems ^ which may be solved (amoiig several 
possible strategies [^) by going over to the so-called second order theories (SOTs) 181. The 


performance of SOTs with respect to RHICs is analyzed in 


19|. 


There is not a single SOT framework as compelling as the Navier-Stokes equations in 


the non-relativistic regime [18, l2Q|, l2l| . However, in the linearized regime they all agree in 


providing a set of Maxwell - Cattaneo equations 


22| for the viscous stresses and conduction 


currents, while they differ in the way the transport coefficients in these equations are linked 


to the underlying kinetic theory description 


23 


29| . For this reason in this paper we shall 


consider only the linearized regime. This is what makes our analysis preliminary, because it 
is likely that the most important effects of the ffuid - held interaction will be connected to 

M], field - turbulence interactions js] and 


nonlinear phenomena such as inverse cascades 


30 


hydrodynamic instabilities 32|. However, as we shall see, already in the linearized regime 
there are significant qualitative differences between SOTs and FOTs, and between SOTs on 
fiat and expanding backgrounds. 


There is a large literature on cosmological models based on SOTs j33l-l36l| . This literature 


focused for the largest part on homogeneous models, where the interest was in how viscous 
effects modified the cosmic expansion and contributed to entropy generation. These analysis 
showed that there are meaningful differences between ideal, first and second order theories 
even at the largest scales. To our knowledge, the application of SOTs to inhomogeneous 


models is less developed than FOTs 


37 


This consideration also contributed to make an 


analysis such as this paper a necessary first step. We note that a family of exact solutions 


for the Boltzmann equation in expanding backgrounds with a we 
limit is known, which provides a helpful test bench for the theory 


defined hydrodynamic 
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In summary, we shall adopt the so-called divergence type theory supplemented by the 
entropy production variational principle (EPVP) as a representative SOT, but will regard 
the transport coefficients as free parameters, rather than attempting to derive them from 


40| . For this reason, our analysis is relevant to any SOT 


an underlying kinetic description 
model. 

The equations of the model are the conservation laws for energy - momentum and charge, 
the Maxwell equations, and the Maxwell - Cattaneo equations providing closure; for a de¬ 
tailed derivation see |^. In the linearized regime, these equations decouple in three sets 
of modes, sound waves, incompressible shear waves, and electromagnetic waves coupled to 
conduction currents. We shall consider only the latter. 

We shall model the Universe during reheating as a spatially flat Friedman - Robertson - 
Walker (FRW) model, whose metric in conformal time is {rj) [—-|- dx^ -|- dy'^ + dz ^], 

a [r]) being the conformal factor. We shall assume for the fluid an equation of state 
p = (1/3) p and vanishing bulk viscosity. Under these prescriptions, FOTs lead to con¬ 
formally invariant equations js]. Therefore the electromagnetic helds are suppressed by 
a factor, on top of the hydrodynamic evolution. We shall consider the evolution of 
electromagnetic helds in an environment where the temperature is higher than the QCD 
phase transition temperature, i.e. a scenario where SOTs seem to correctly describe the 
state of the matter. 

Unlike FOTs, the equations derived from SOTs are not conformally invariant: the ex¬ 
pansion of the Universe creates temperature gradients which couple to the huid velocity and 
conduction currents. This leads to a weaker suppression of the magnetic helds than expected 
from a FOT framework. This is the main conclusion of this paper. The ehect is not large, 
but suggests that these SOTs models may be more sensitive to nonlinear ehects, such as 
hydrodynamic instabilities, than FOTs or even SOTs on non-expanding backgrounds. This 
possibility will be investigated elsewhere. 

The paper is organized as follows: In Section [TTl we introduce the formalism and the 
covariant equations of second order hydrodynamics. We analyze the conformal invariance 
of the theory and derive the equations for the helds as well as for the viscous stress and 
conduction current, showing that the latter are explicitly non conformally invariant. In 
Section IIIII we linearize the equations and propose a simple, toy model, to solve them. In 
Section HV] we consider the homogeneous case /c = 0, that permits to study the electric held 
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and conduction current separately from the magnetic field. In section |V] we consider super¬ 
horizon modes of astrophysical interest, i.e., /c <C 1, and hnd that the magnetic helds evolves 
in a way clearly different than in FOT’s models. In Section |Vl] we summarize and discuss 
our results. We leave for the Appendix the analysis of sub-horizon modes fc S> 1 as they 
are not as astrophysically interesting as super-horizon modes. In Appendices [B] and O we 
quote some secondary results and technicalities for the reader interested in those details. We 
work with signature (—, -|-, -|-, -f) and natural units h = c = ks = thus time and length 
have dimensions of energy~^, while wavenumbers, mass and temperature units are those of 
energy. 


II. GENERAL RELATIVISTIC FLUID EQUATIONS 


A. The equations in covariant form and their 3 -|- 1 decomposition 


We consider a system composed by a neutral plasma plus electromagnetic held in a hat 
FRW universe, whose metric in conformal time is ds^ = [r]) [—47]“^ -|- dx"^ + dy'^ + dz^], 

a {rj) being the conformal factor. This form of the metric is obtained from the one written in 
physical time t by dehning dr] = Hodt/a{t), with Hq the Hubble constant during Inhatioql. 
If for Inhation we consider the de Sitter prescription, ai{t) = exp (Hot), then aj (g) = 
1/ {1 — g) with g < 0. If for reheating we accept that during that period the Universe evolves 
as if it were dominated by matter |ll|, then aji{t) = (1 -|- (3/2)and consequently 
Or (h) = (1 + ^/2)^- Observe that we have matched the two expressions ai t = g = 0 such 
that 0 /( 0 ) = a/j(0) = 1. As Hq is a hxed, characteristic energy scale, we can use it to build 
non-dimensional quantities, as we did with conformal time, e. g. we dehne dimensionless 
lengths and corresponding wavenumbers as I = HqE and k = k/Hq. Magnetic and electric 
held units are energy‘s so we write B = B/H^ and E = S/Hq. To complete, we quote the 
temperature T = T/Hq and the electric conductivity = ad Hq. We use greek letters 
to denote space-time indices, and latin letters when we deal with spatial-only components. 
Besides, we use semicolons to express covariant derivatives and commas to denote partial 
derivatives; in particular a ’prime’ will denote partial derivative with respect to conformal 
time, i.e.. A' = dA/dg To evaluate the diherent covariant derivatives we need the Christohel 


^ With this definition, y is already dimensionless. 
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symbols, whose only non-null components are Fqq = a'/a, F?- = a'/aSij and Fq^ = o! jaS'y 
Let be the fluid four-velocity. We decompose it as 

= 7 (f/^^ + ( 1 ) 

with 7 = \J\ — v^. It is satisfied that = — 1 and = 0. is the velocity of 

fiducial observers and represents deviations from Hubble flow, i.e. peculiar velocities. 
Each of these velocities defines a congruence of time-lines, for which there is an orthonormal 
space-like surface defined through the projectors 

The matter is described by the energy momentum tensor, which we decompose as 

with 

= (^p + p)uV + pg^’' ( 4 ) 

and 

t"" = TtJ-iC"" (5) 

15 

the viscous stress tensor. In eq. (|5]), r is a characteristic relaxation time and is a 
Lagrange multiplier whose evolution equation will be given below; for r —)■ 0 it reduces to 
the FOT dissipative shear viscous tensor. We write the electromagnetic held tensor in 
3-1-1 form relative to the fiducial observers as 

— Ay^^ = U^Ey “ Ey E TjfiydpU^B'^ ( 6 ) 

with 7 O 123 _ \jl^t {—g^y)]~^^‘^. For future use, we define Observe that the 

electric and magnetic fields are obtained from ([ 6 ]) as = F^’^Uy and = {l/2)g^''°‘^UyFap 
respectively. The electric current is 


.P = (7) 

with 

T" = JTT 2 C (8) 
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where is another Lagrange mnltiplier whose evolntion eqnation is also given below, and 
that for r —)■ 0 gives the nsnal Ohm’s law. 

Althongh we shall regard J^n in eqs. (ED and ([HD as free parameters, we observe that these 


eqnations may be derived from a linearized Boltzmann eqnation 
seen to be 


4l| . in which case they are 


^n 


Dpfl 


\ 1 ^ 
-U Px\ 


(9) 


with /o the one particle distribntion fnnction. Dp = (27r) 2d'^p5 {p^Pix — m?) the integra¬ 

tion measnre (m is the mass of the plasma particles), and where F is a mnltiplicative factor 


or F are Marie’s prescription 


42 


m 


4^,1^ whereby F = \u^p^\. 


in the linearized collision integral. Common choices 
i.e. F = const., and the Anderson-Witting proposal 

Observe that in eq. ([6D we dehned the electromagnetic field relative to hdncial observers. 
It is also with respect to this velocity that we shall dehne the ‘total time derivative’ or ‘dot 
derivative’, namely . The ‘total spatial derivative’ is accordingly dehned as 

The eqnations we have to solve are the conservation eqnations (matter conpled to the 
electromagnetic held pins charge conservation). Maxwell eqnations and two eqnations that 
describe the evolntion of the Lagrange mnltipliers and The conservation laws are 




( 10 ) 




( 11 ) 


and Maxwell eqnations in covariant form read 


F^\, = -Jf^ 


( 12 ) 


= 0 


(13) 


To onr pnrposes the best is to rewrite the previons eqnations in 3-1-1 form relative to 
hdncial observers. This is achieved by projecting each set along and onto its orthogonal 
snrface described by The projection along is dehned as [4^ — 

and the one onto the orthogonal snrface as For the set ffTop we hrst replace 
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expression ([T]) in eq. (jl]) and define 


P = 7 ^ (P + P) - P 
P = = P + 3 ( 7 ^ - 1) (P + P) 

= 7^ (p + p) v'" 


= ^2 ^ti^u _ f ^2 _ 


We thus write eq. ([3]) as 


15 


For eqs. ([7]) plus (jH]) we directly obtain 


= p,7 {U^ + v>^) + -tT^C^ 
o 


(14) 

(15) 

(16) 

(17) 

(18) 

(19) 


To find the evolution equation for the plasma we assume the equation of state p = p/3. For 
the projection along of eq. (HB we have 


1 [(47 - 1) + i P'W).,. + ~ [(47'^ - 1) p] + ^ 


= E’' ( JPqV^ + 


while for the spatial projection we obtain 


( 20 ) 


p,„A*- + 1 Yi^pvV] AS + 57 ^"a; + i^'AS 


3 " ^ - Q 

-^r'ry“"AS - jTT*T„C“'AS - 7rTyyAS 


p^E'^ + RP p^v^ + tTC 


a 


( 21 ) 


For eq. (HU using (IT^ we have 

'yPq,flU^ + lPq,^pV^^ + Pq [ 7 ^^^“^;^“^ + 7^^^ + 

+eVTy7'‘ + "-tT^Z = 0 


( 22 ) 


As Maxwell equations are already written in terms of the projection is straightforward. 
For the inhomogeneous Maxwell equations ffT2|) we have 


q: = p,- 

ASB” = - 24 a;B'' + AS7 ”‘’'”t;,B„„ - AS J“ 


(23) 

(24) 




while for the homogeneous ones ffT^ we obtain 


B-p = 0 


1 . „ , a' 


+ 2—a^:b^ + A^:w = 0 


2 7 


(25) 

(26) 


We now discuss the equations for the Lagrange Multipliers and see j40| and [41 1 


for details. C ^md Cy are transverse with respect to and Ay is also traceless, i.e. they 


satisfy 


CX = o = C-w^ C = o 


(27) 


Their evolution equations in covariant form are s 


spending Minkowski expressions given in Ref. Mj. We obtain: 


xaightforwardly obtained from the corre- 




(28) 


Ai [Cf,u + rh“h^Ca/3;7«^] = 


A4. 

A5 T 

7 T 


tA^C 

■^5 rjn2 S) 






(29) 


with = (1/2) [u^-u + Ui,-^] — {l/3)u°^h^i, the shear tensor. In the derivation from linearized 
kinetic theory the functions An are given by 4l| : 


An= L>p|-Map"r/o 


(30) 


We only mention this because it makes it easy to check the dimensions of Bn and An] oth¬ 
erwise we shall regard them as free parameters. The dimensions of the different expressions 
under the integrals are [/o] = 1, [Dp] = [uaP°‘] = E, with E meaning ’energy’ and 

consequently [An\ = and [Bn\ = E'^^^. As the only energy scale of the plasma is its 

temperature, we rewrite eq. (|1]) as 

(31) 

15 

and eq. ([8]) as 

= ^C2tT^A (32) 

o 

with Cl, C 2 dimensionless, 0{1) coefficients. 
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B. Conformal Invariance 


To analyze conformal invariance we begin by rewriting the coefficients in eq. fl28l) and 
as 

A 2 _ hi ^i _ 1 1 _ &3 _ J rj, /oQ\ 

^3 “ T’ ^3 “ A “ T 3 ’ A “ ' 

where bi, 62 , bs, di are again numerical, ~ 0{1) coefficients. Therefore the mentioned eqs. 
read 


M 




T 

,/3 


e2T3 


b/3'V 




diT^^^Cliv = [C^lV + 


di 


u 


+ Cau(^'^ - 


We now transform the different quantities in the model according to 


(34) 


(35) 


M'' 


= - - )■ u^ = aUf, h'""' = — -^ -)■ = aa^ 

CL Oj 


flU 


(36) 


(and similar rules for and A^‘' 


C" = ^ ^ C, = C"" = ^ ^ 


(37) 


and 


pfiu 

p = p ^ p^^’^ = _ 


P P Pq ra rr To 

d=-7, P=-7, Pg = -7, J =-X^ T = —, T = aT 


^ 4 ’ ry 4’ 4 

a 


Replacing these transformations in eqs. (fT8|) . ([7]), fl3lD and 


we hnd 


rpfllP _ -^0 /iiy _ _ T/i _ ^ 

“ 0® ’ a6 ’ a* 


(38) 

(39) 

(40) 


and the set of eqs. (l20|) - fl2T]) becomes 

^ ( 47 ^ - 1 ) P' + + 7P7j^^ + + 5fT^Tj<°^ + rT^Cj 


=h;^ 


7PgWj + J^T^Cj 


(41) 




P,E^ + t\B^ { P,V, + jff% 


(42) 
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while for eq. fl2^ we have 


IPq + + Pql 




Q:,/i 


m 


+ Pqiv% + e^fT^T^ 


+ -fT^C^, = 0 


(43) 


It is a well known result that Maxwell equations are conformally invariant. For the homo¬ 
geneous equations it is a trivial result, and for the inhomogeneous equations it is directly 
apparent from the transformation law for and the last of exprs. (1401) . Therefore trans¬ 


forming eqs. fl2^ - fl2B]) we have 

y. = p,+ df fy» (44) 

E“ = - p,v‘ - jff^C (45) 

P, = 0 (46) 

B" = -e‘\Ej (47) 


Notwithstanding, when we apply the above conformal transformations to the evolution 
equations for and conformal invariance is lost. To see this, we replace and 
from eqs. ([1]) and ([6]), and use the conformal transformations defined above to obtain 


4 = 2^7 


r ^ ^ ^ ^ 1 


a' ~ _ o ~ 

Efj_ -|- F 

- 

Ca,0 Co T n Co,(3 

a 


(48) 


and 







(49) 


-r— 7 ' + 7 3- + n(; 
7 V a ’ 


c -fP 




In both equations, the terms proportional to a'/a do not cancel out and this fact makes 
the two equations non conformal invariant. As the fields evolve coupled to this plasma, the 
conservation of the magnetic flux during their early evolution is lost. To have a glimpse of 
the effect of this coupling on the amplitude of the magnetic field, we shall solve the equations 
in the linear regime. 


III. LINEAR EVOLUTION 

The system of equations that describe the evolution of the plasma is non linear. We shall 
study the linear regime, that is suitable for small amplitudes. We shall also consider that the 
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plasma is neutral, i.e., we assume pq = 6pq = 0. First order quantities are , 6p 

and the electromagnetic field. Writing H (p) = a'/a, the linear equations read 


Sp' = -gPohj- 

^ -i/5,. + 


V, = 


c,= 

C' = 


4po 

H(V)- 


10 Po 
4 1 

( T 


c 




b2T 


Cl + , ~rp Ei 
02rTn 


B^' = 


jfT^C 


k ,j 


(50) 

(51) 

(52) 

(53) 

(54) 
(56) 


where we see that at this level the plasma equations have separated from the electromagnetic 
equations, so from now on we concentrate only in the latter as our focus is the electromagnetic 
field evolution. Before going on, observe that if we set r —)■ 0 in eq. fl53l) we have that Q = 
(6i/To) Ei. Replacing this expression into Ampere law, eq. (|54|l . the factor that multiplies C 
in the last term of the r.h.s. becomes (e^/3)fTQ feii?*, and recalling the constitutive relation 
between electric field and density current, J* = dcE^, we can read the expression for the 
(commoving) electric conductivity: 

(56) 

Observe also, that due to the conformal scalings fl3^ the physical and commoving electric 
conductivities are related in the usual way, i.e. dc = d'c/a. We then rewrite eq. fl54|) as 

(57) 

To go on we change the time dependence from p to u = (1 + p/2), whence d/dp = 
{l/2)d/du and H = 1/m. Assuming incompressible evolution and transforming Fourier we 
get 

dC (k, u) 
du 

= 2ie"\yB'^ (fc, u) - (^, u) (59) 

dR* (fc, u) ■■ ■ -i. N 

—^ = -2fe«,F£)*(A/,q (60) 


1 1 

62 ^ u 


Oh 

C (k,u)+^K (ku) 


(58) 
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We shall not attempt to solve system fl58|) - (l59|) numerically, as this would oblige us to 
stick to a specihc range of parameters. Instead, to have a glimpse of how the system behaves 
we assume a simple conhguration given by 


k = kz, Bi = Byij, Ei = E^x, 0 = CxX 


(61) 


Dehning the matrices 


A = 


Ex 

\By J 


and 



/ 1 

bi 1 

0 ' 


/ 1 

0 

0 ' 


b2T 

62 xTq 


U 

- 

o'cTo 

bi 

0 

ik 

, H = 

0 

0 

0 


0 

ik 

0 


0 

0 

0 


(62) 


(63) 


the system of equations for the electromagnetic sector can be written in matrix form as 


A' + 2HA = 2HA 


(64) 


where now a ’prime’ denotes derivative with respect to u, i.e., / = d/du. In spite of its 
simple form, it is rather difficult to solve eq. flMl) exactly, except for the homogeneous 
mode, A; = 0. We begin by solving this case and then consider perturbatively the case 
A: 1, that corresponds to modes well outside the particle horizon as e.g., the galactic 

scale. The solution for modes A: 3> 1 is given in Appendix lAl 

To appreciate the features of the SOT evolution, it is convenient to keep in mind their 
behavior in the r —>■ 0 limit, whereby the model reduces to a FOT. In that case system 
(I58|l - (j6nil plus (j56il and model (IbTll reduces to 

E {k, u) 


du 

B^{k,u] 

du 


= 2ikB — 2acE 
= —2ikE 


(65) 

( 66 ) 


and this (conformally invariant) system can be combined to give a wave equation whose 
solutions are the exponentials with 7 (±) = dc/2 ± — k'^. Observe that when 

A; —)■ 0, 7 (+) —)■ (Tc and 7 (_) —)■ 0. The second solution describes the “frozen” magnetic held, 
and the hrst the “discharge” of the electric held due to the resistivity of the plasma. If A: 7 ^ 0 
we have the well known pure exponential decay. 
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IV. ANALYTIC SOLUTION FOR THE HOMOGENEOUS MODE R = 0 


In the fc = 0 case, eqs. fIMl) may be solved in closed form. We then begin by pntting 
/c = 0 in matrix H and the r.h.s. of eq. flM|l eqnal to zero. Proposing as solution a time 
dependence of the form A* (t]) = A^°^*exp (—and imposing that the determinant of 
the resulting system be zero we obtain the eigenvalue equation: 

1 


;^(0)2 


62 r 


_A(o) 


= 0 
02 'r 


whose solutions are 


\(o) _ n 

"^(0) “ u 


x(o) _ 
^(±) “ 


2&2'r 


(^1 ± - 462 ro-cj 


(67) 

( 68 ) 

(69) 


Observe that there exists a critical relaxation time, fc = 1/ (4625‘c). Also, and more impor¬ 
tantly, when f —)■ 0 we have that A(_) —)■ dc while A(+) blows out. Therefore A|q| and 
converge to the roots of the FOT model. The corresponding eigenvectors are 


A(o) - 


/o^ 

0 

VV 


a(o) - 


( 1 \ 

o'cTp 


blA(+) 

° / 


a(») - 


1 t 

O'cTq 


6iA(_) 

0 


(70) 


” / 


To hnd the solution of the inhomogeneous equation we propose 


A* = 


(m) aJJI + a(+) (n) Aj^^.e + «(_) (n) aJ^^.c 


0 ( 0 ) 


ho) 


h+) 


(-)' 




(71) 


and substitute in eq. fl6T)) . For a(o) {rj) it is straightforwardly obtained that a(o) {rj) = const. 
Recalling that this coefficient corresponds to and that this eigenvector represents the 
magnetic field, this means the obvious result that the commoving field remains constant and 
consequently the physical magnetic intensity will decay as cx [t]). The other coefficients 
satisfy 


®(+) - 


a(-) - 


A^°^ 9 r 

^(+) 2 r 

AA(o) u 

A^°^ 9 
A(-) 2 


Cl(+) -|- 0(—)C 
_ ^-2AA('')« 

a(+)e 


2AA(0)n 


AAWn 

with AA^°^ = A®^ — System fl72|) - fl73l) can be reduced to 


'(-)• 


«(+) = 


AA(°) 


u 


(0) 2 


lV) 


i(-) + a(_) 


2AA(“)« 


(72) 

(73) 


(74) 
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plus an equation for a(_): 




+ 


2AA(°^ - 


u 


H-) 


x{o) 

+ 4^a(_)=0 
u 


(75) 


which through the change of variable 2 ; = —can be rewritten as a Kummer equa¬ 
tion, whose solutions are the Confluent Hypergeometric functions 47|, |48|]. Two linearly 


independent solutions of this equation are 


m 


^ = U (^2A|°yAA(°), -1, -2 AA(°)m) 

and ^—1 — 2A|°^^/AA*^°\ —1, 2AA^°^nj. For the given parameters, both 

functions converge for u —>■ 0 and u ^ 00 ^|. Therefore we write a(_)(M) = aa|^yM) -f- 
/9a(_)(w), with a and (3 constants to be determined by the initial conditions. 

Before analyzing asymptotic behaviors it is important to establish the (conformal) time 
interval where the evolution takes place. As said before, we are considering conduction 
currents, which are likely to be made of the lightest supersymmetry particle s-t. This 
means that we are considering times before the establishing of the standard electron-proton 
plasma, which we can estimate as being the time of the QCD phase transition. Moreover, we 
are interested in the final states of the magnetic evolution coupled to this current, because 
it would give the initial conditions for the subsequent evolution of the held in the standard 
proton-electron plasma. If we take the standard value of the Hubble constant during Inhation 
H = 10^^ GeV and the Planck mass as rupi ~ 10^® GeV, we estimate the temperature at 
the onset of reheating as T^h = ~ 10^® GeV. This plasma cools down due to the 

expansion as = u~‘^. The electroweak phase transition took place at a temperature 
scale of Tew ~ 10^ GeV, therefore the dimensionless, conformal time elapsed since the 
onset of reheating can be estimated as Aui^r^ew) — \/Trh/Tew — 10® S> 1. Moreover, if 
we consider the QGD phase transition, for which Tqcd 10 ^ GeV, then Au(r_qcd) - 
\/Trh/Tqcd — 10 ® 3> 1 . Therefore to hnd the sought initial values for the subsequent 
evolution in the radiation era, we can safely take the limit u S> 1 ~ cx) throughout. 

When n —)■ cx) the Gonfluent Hypergeometric functions can be always approximated as 
471. 1481] f/(a, 6, z) ~ and as e e we get 


C (m —)■ 00) ~ a (— 2 AA^®^m) ^ 


-2A™/AA(0) 




(76) 


and 


E (u ^ 00) ~ a 


T 


0 CTc 
bi A(_) 




(-)‘ 


(77) 
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We see that ( oc (61/To) E. Observe that for f > fc the solution becomes oscillatory. This 
behavior has no analog in FOT’s. However, if we take the limit r —)■ 0 in eq. ( 1771 ) . where 
A^^/AA® —)■ 0 and dc, we recover the FOT results. 

V. NON-HOMOGENEOUS MODE K < 1 

Although the homogeneous mode has the appeal of affording a full analytical solution, it 
is clearly not very interesting from the cosmological point of view. In this section we shall 
consider the set of modes which are most relevant to cosmology, namely modes which are 
far beyond the horizon at reheating, k 1. Now the magnetic held is no longer decoupled 
from the electric held, and we expect to hnd some feedback from the latter on the former, 
eventually reducing the cosmological suppression and the exponential decay found in 
the FOT analysis. We shall not attempt a full solution, but rather analyze the asymptotic 
behavior of the magnetic held. 

To solve for A; 7 ^ 0 we begin by solving perturbatively the eigenvalue equation det H = 0. 



(78) 


by proposing 


A(‘) = A'“) + + ■ ■ ■ 


(79) 


After replacing (17^ into fl75]l and keeping terms up to order we hnd 



(80) 


a(0) 

^( 2 ) _ 



(81) 


with A|°^^A|°^^ = <Jc/{h 2 T). To hnd the eigenvectors we again set to zero the r.h.s. of eq fIMl) . 
and propose the new eigenvectors as linear combinations of the /c = 0 ones, i.e.. 


A<''> = o,(,)(A^)AjJ’ + 0(+)(A^)Aj“’, + o,_)(A^)Aj“’, 


(82) 


The results are shown in Appendix |Bl 

To solve the time evolution we rewrite eq. fl6T)) as 



du 


(83) 
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In the above equation, upper index in S denotes row while lower index denotes column. 
Keeping in mind that the physical range of u starts at k = 1, we Laplace transform A* as 


F\s) = / du 


and so eq. 


becomes 


d 


— [sF\s) - A'(0) + 2EiF^{s)] = 2SidjF{sy 


As the term involving the initial condition vanishes upon deriving we obtain 

d 


ds 


[e(s);F^P)] = 


where we have dehned 

Q{sy = 2^ + sS} 

We now introduce the inverse matrix of 0(s)*, M{s)y i.e. 

Q{s))M{s)i = SI 

and dehne a new variable such that 

F^s) = M{s)iK{sy 


(84) 

(85) 

( 86 ) 

(87) 

( 88 ) 

(89) 


Replacing in eq. (j86ll we obtain the following equation for iL(s)*: 

A [A'-(s)] = -26lM(s)lK(s)'‘ (90) 

We see that for i = 2, 3 the solutions are constants. For i = 1 we have 

-+ 2M}(s)K\s) = -2 MUs)K^ - 2MUs)K^ (91) 

ds 

Previously, we have solved the eigenvector equation for the homogeneous system, i.e., 
with a = 0, +, — (no sum over greek indices). To avoid cumbersome 
notation from now on the label {k) is omitted. There exists the inverse matrix to A(^^, 

i.e. 

nf(92) 

that also satishes 

= s‘ (93) 
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Using this result we can write S*- as 

E-=5^A,„,A;„,nf' (94) 

(a) 

from where we can write the matrix Mj as 

U(») = !]''(«) + 2^1"))nf’ (98) 

Q; 

We now solve eq. fl9Tl) . The homogeneous solution is straightforwardly obtained and reads 

Kk^) = n (96) 

a 

with 

A,„, = (97) 

(no sum over a). Observe that using relation ([93]), A(^a) satishes 

= l (98) 

(a) 

To hnd the inhomogeneous solution we propose K}{s) = L{s)Kl^^{s) and after substi¬ 
tuting in eq. (I9T|) we hnd the following evolution equation for L{s): 

^^L{s) = -2 [M^{s)K^ + Ml{s)K^] J] (s + (99) 

a 

Up to here, all the developments have been exact. However, to hnd solutions that rep¬ 
resent the resulting held after the evolution in the reheating plasma, we must solve (199]) in 
the asymptotic range s —)■ 0 (i.e., u —)• cx)). 

A. Solutions for s ^ 0 (u —)■ oo): 

We now look into the small s limit. We begin by recalling that one of the eigenvalues, 
A(o) goes to zero as A; —?• 0, while the other two A(+) and A(_) remain hnite. Therefore for 
small enough s we can take s A(±), but we cannot assume s < A(o). Therefore, retaining 
this last eigenvalue explicitly, we get (up to an unessential constant) 

/qo™(s)« (100) 
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and 


A^ (s); = ([2S]-‘)1 + Ajn” (s + 2A,„,) - (2A,„|) 


( 101 ) 


with 


( 62 T 


[ 2 : 


1-1 


^ 0 

2 ^ 


ibi 


‘IkT'o 


\ 


0 

ib2TTod-c 
y 2bik 


0 —- 

^ 2k 


( 102 ) 


'2k 2fc2 / 

The solution of the inhomogeneous equation now reads 

-2 


L{s) = 


2An + 1 


(|2H]-‘),-A;„,nf (2A,„,)-‘ 




([2Sr‘)‘-A‘ nf (2A,„,) 


-1 


— [AlU^K^ + AlUlK^] {s + 2Xo 


2Ao 


I (s + 2 Ao ) 
+ K^ 


2 A 0+1 


(103) 


whereby 

{s) = L{s)Khom {s) 
-2 


2^40 + 1 


Al Tt(0) 

([2=]-): - 


2A 


(0) 


A 


( 0 ) 


A;„,nf + AjojnfA’ 


+ (s + 2 Ao ) 


Al 


3 

-2A 


2A 


( 0 ) 


( 0 ) 


1 (s + 2A(i))) 
(104) 


The different functions then read 


F‘ = M}(s)kHs) = (pS]-'), A-‘(s) + ([2S]-‘)^ A-2 (s) + (|2H]-‘)3 A'=>(s) 

+ i - Aj, [n;A‘ + n»A'= + n»A’] (106) 

_s + 2A(o) 2 A(o)J 

Our main interest is as it is directly related to the magnetic field. The calculations are 
long but straightforward and are shown in Appendix O The result is 

^on? 1 r . 1 T-r{0) , A 1 Tt{ 0) r^sl , ^on? 


F ~ 


■7-Try;— Afoin® AT^ + Al.uf^K^ 

(s + 2Ao) A(o) L ' ( 0 ) 3 

^0^3 ^3 


+ 


> + 2Ao)^+^^‘°) 




(106) 


(s + 2 Ao ) (s + 2 Ao ) 

The calculation of elements If® is also rather long but straightforward, here we quote the 
one in the term with as this term gives the main contribution. It reads fl)* = ikb 2 fTQ/bi, 
and we then have 

b2fToK^ 


F ~ 


6i(s + 2Ao)^+'^(°) 


ik 


(107) 
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where A(o) — b 2 Tk‘^/ac > 0. The corresponding anti-transformed function is 




h2rToK^ .. _ /1 


hi 


u 


exp [—2 Aom] 


( 108 ) 


Observe that due to the presence of the factor the magnetic held decays slower than 

the exponential law of FOTs, even at large times. 


VI. CONCLUSIONS 


In this paper we have studied the evolution of electromagnetic helds coupled to conduction 
currents during the reheating era, using second order causal hydrodynamics to describe 
the evolution of the currents. The evolution of the magnetic held occurs well before the 
EW phase transition, during an epoch where the standard proton-electron plasma is not 
established yet; the conduction currents we consider are likely to be made of the lightest 
supersymmetric partner s-r. The main motivation behind the choice of SOTs is the well 
known fact that hrst order theories (as e.g. relativistic Navier Stokes equation) have severe 
problems of causality and have no stable equilibrium states. Also, SOTs behave quite well 


at describing RHICs 


13 


19| . where a plasma much like the one in the very early Universe 


is supposed to be created. Thus, although there is not a preferred SOT framework yet, 
it is important to begin to study different plasma effects in the early Universe using those 
formalisms extended to general relativity. We adopted the so-called divergence type theory 
plus the entropy production variational principle (EPVP) as a representative SOT, but 
regarded the transport coefficients as free parameters, rather than attempting to derive 


them from an underlying kinetic description 


40| . In this sense our analysis is relevant to 


any SOT model. When extended to General Relativity, we found that the resulting theory 
is not conformally invariant: the Maxwell-Cattaneo equations that describe the viscous 
stresses and conduction currents lost this symmetry. As these equations are coupled to 
Maxwell equations, the consequence is that the magnetic flux is not suppressed by the 
expansion as quickly as in the Navier-Stokes theory. This might provide higher intensities 
as initial conditions for the subsequent evolution during radiation dominance. The physical 
explanation is that the expansion of the Universe produces temperature gradients which 
couple to the current and generally oppose dissipation. 

To pursue the analysis we considered only the linear evolution because in this regime 


20 














all SOTs agree in providing the set of Maxwell - Cattaneo equations. Our goal was to 
identify the qualitative differences between FOTs and SOTs, this is the reason why we did 
not attempt to give numerical estimates of the resulting amplitudes. We have found that 
the held decay in the homogeneous mode may be oscillatory. Even in the purely decaying 
regime, for inhomogeneous modes there is a power-like correction to exponential decay, with 
a positive exponent. This suggests that the unfolding of hydrodynamic instabilities in these 
models follows a different pattern than in hrst order theories, and even than in second 
order theories on non expanding backgrounds. The study of the non-linear hydrodynamic 
instabilities is the next step in the research of primordial magnetic helds evolution within 
SOTs. 
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Appendix A: Large k modes 


Although small scales are of little astrophysical interest concerning galactic magnetism, 
for completion we devote this appendix to analyze their evolution with the formalism con¬ 
sidered in the paper. Moreover, in this case the mathematics is much simpler. We shall see 
that the effect of the relaxation time r is to add damping, while the effects of conformal 
invariance breaking is to add a slight amplihcation of the magnetic held if the temperature 
(and therefore the conductivity) is low enough. The equations were 


dC 


2 

&2T 


2 

u 



2bi 

hfTo 


E 


dE 

du 


+ 2ikB + 


2d cTq ~ 


——h 2ikE 
du 


0 

0 

0 


(Al) 

(A2) 

(A3) 
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It is convenient to introduce £ = E + B and B = B — E to get 


Now we write 


du 


2 

&2'r 


c {£ B) = 0 

2b2fTo ^ ^ 

(A4) 

d£ , ^ 2dcTo ~ 

^ +2ik£+ ; °C = 0 

du bi 

(AS) 

f mB = o 

du bi 

(A6) 


S = 

B = BoE^^^ 
( = 


(A7) 

(A 8 ) 

(A9) 


with the understanding that the pre-exponentials are all slowly varying functions of time. 
Collecting positive and negative frequency oscillations we get 




2 2 ' 

2ik -|- - - 

02'r u 


C+ + r^Bo = 0 


du 


&2'rTo 

25'cTo 


C+ = o 


and 


du 


, 2 2 

-2ik -|- -- 

O 2 T u 


C- 


bi 


&2rTn 


-^n = 0 


Leading to 


d^Bo 

du"^ 

d^So 


+ 


2ik + 

02 T 


du^ 


+ 


-2ik -f 


62 7 


d£o 2(TcTo ~ 

+ —^C- = 0 

du bi 


dBo ac „ 2 dBo 

-r + rBo = -^ 

du 02 u du 

d£o ^ ^ d,£o 

du 62 ° M du 


(AlO) 

(All) 

(A12) 

(A13) 

(AM) 

(A16) 


Let us analyze the equation for Bo- Setting the r.h.s. of eq. (lA14p to zero, the solutions are 


with 


The roots are 


+ 


2i' 
2k- — 
02 T 


w- ^ = 0 
02 T 




±1 




(A16) 


(A17) 
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The slowly varying solution being uj+. Therefore we postulate 


i3o = 


to get 


da 

du 




+ = 0 


du 

which for b a becomes 


du 

^ (a;+ae*"'+“ + 

(Jj Uj 


U 


da 


2a;_i 


du (a;+ — U-) u 
db —2u 




du (a;+ — U-) u 


whose solution for a is 


a = u'^ with a = 


2ca_i 


Therefore we get 

Bq ~ exp iu+ 

When k is very large we have 

a;+ = 


-2i 


la, 


a;+ — U- 

1 


a;+ — a;_ 


In M + M 


so 


and 


therefore 


We may write 


where 


Re 


Re iu+ = 

2u+ 

— a;_ 

2ca_L 


2 1 + ikb2f 

-dc 1 
2 bnf^k"^ + 1 


b2Tdc 


[1 + ikb2ff 

blr^k^ - 1 


= b2fdc 


u+ — U- + 1]^ 

\Bq\ ~ exp {A [ticlnti — u]} 
blf^k^ - 1 


Uc = 262T 


bnf‘^k‘^ + 1 


Z\ = 


(Jr. 


2 b%T‘^k?‘ + 1 


(A18) 

(A19) 

(A20) 

(A21) 

(A22) 

(A23) 

(A24) 

(A25) 

(A26) 

(A27) 

(A28) 

(A29) 

(A30) 

(A31) 
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Bo grows up to Uc with an amplification factor 


Bo (uc) 


= exp {Auc [In {uc/ui) - 1 + Ui/uc]} 


Bo {ui) 

Of course, provided Uc > Ui. For all practical purposes, the amplihcation is 

I Bo (uc) 


(A32) 


Bo (ui) 


1 dc 


(A33) 


Appendix B: Eigenvectors for /c <C 1 

After long but straightforward calculations the eigenvectors for the perturbatively cor¬ 
rected eigenvalues are: 


/ b\ik \ 


^(0) - 


CTcTo 

_ ih. 
1 


\(k) _ ^ 
\o) - ~ 


Or. 


\ ^ ) 


A(fe) 


1 - 


CTeTb 


biX 


(0) 

(+) 


aa(o)2 

^(+) 

AA(0)2 

ik 


O^cTq 


bi\ 


(0)2 

(+) 


a(0 _ 


1 - 


ucTp 


biX 


w 

(-) 


x(0) 

^(+) k^ 

.(0) aa(0)2 
^(-) 

-1 

^ “ AA(0)2 




AAh—’i Jc 

_ oth 


feiA 


( 0)2 

{-) 


\(^) _ \(0) I 

A(+) - + 




,( 0)2 ,( 0 ) ,( 0 ) 

^(+) \+)\-) 


x(0 _ x(0) , 




,( 0)2 


x(0) ,(0) 


(Bl) 


(B2) 


(B3) 


Appendix C: Solving for 


Explicitly we have 


F'^ = 
+ 


([ 2 : 

([ 2 : 

([ 2 : 


.-lo Agn? 

2Ao 

,_1,3 A3n° 
2Ao 

A^no 


D - 

J /3 


t1\3 


) - 
/2 


2Ar 


s F 2Ao 

A^n° 

^3 ^ 

s -|- 2Ao 

2 Amo 

^ 

s F 2Ao 


(Cl) 
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Replacing the different expressions we obtain 


2Ao |2Alo + l 


1 tt(0) 


Al tt(0) 


+ (psi-);- 


Al n'^ 

-1x1 


(s + 2A(o)) 


A;„nf a' 2 + A|(„nf A=>] + 


^(0) L ' ' 

Agn; [ -2 [ 

(s + 2Ao) ^ 2^0 + 1 

(Psl-’S-Mi 


{s + 2Ao)"^(°^ 


1 1 Al .nf 
\K^ (s + 2A(o)) 


Al(o) ^ ^ J + (. + 2Aof 

+ [(P=| >3 - A" + +1 


-1\3 Agllg I 3 7V-3 


'3 2Ar 


+ 


Amo 


(s + 2Ao) 


A3+[(|2Sn’-|M|/f3 


A^no 

^^ 0 - 11-2 j ^2 


(s + 2Ao) 

To find the corresponding time dependent function, observe that we can write 

r du (C3) 

Jo {s + ‘^\o)T 

If a —)■ —n the integral diverges for m —)■ 0. But as we are interested in the late behavior 
of the helds it is legitimate to compute the limit when a —)■ —n and discard the divergent 
term (that corresponds to times out of the interval of validity of the approximations made 
in this paragraph). This we do by adding an ’infrared’ cut-off. We then have 

Jn{s)= [ dMM-"+^e-("+2m))“ = r(l-n + e)(s + 2A(o))"”'”' (C4) 

Jo 

where e a small parameter. Developing in Laurent series around the pole we have 

^ ... Ti , /^^ + 2A(0)M [1 , ,, , , , J / , xn-l 


Jnis) ^ 


1 -I- eln 


4-(n-f 1) (s-|-2A(o)) 


s + A(o) 


7- (n -f 1) (s -|- 2A(o)) 


with /i a renormalization constant and ip = T'/T. Finally, for s —)■ 0 we have 


-n+e ^-(s-|-2AfQN)M 


(0)J“ ^ 


7- (n -|- 1) (s -|- 2A(o)) 
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To apply this result to eq. flC2j) we observe that A(o) — b 2 fk‘^/ac -C 1, and therefore can be 
discarded in front of 1. The Laplace anti-transformed different terms that appear in expr. 
flC2p can then be approximated as 


(s -|- 2A(o)) —)■ 


In 


(0) 


/i 


+ V’(i) 


-1 


g-2A(o)ii 


(s + 2A(0)) ^ ^ 


const —)■ — 
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In 

In 


2A 


( 0 ) 


/i 

2A(o) 


-1 




+ V’(3) 


+ ^(2) 


g-2A(o)M 


-1 


g-2A(o)M 


(C7) 

(C8) 

(C9) 


and we see that the contribntion that gives the slower decay comes from correspondence 
flC7ll . We then keep only those terms, obtaining 


F'^ ~ 


Ajn; 


(s -|- 2Ao) ^( 0 ) 


Ajojnf A'" + A;„)nf A'’ 


+ 


A^n? 




-K^ 




-K^ + 






-K^ 


(s -|- 2Ao) (s -l- 2Ao) 

The calculation of elements fl* is rather long but straightforward, here we quote the one in 
the term with as this term gives the main contribntion. It reads II^ = ikb 2 fTQ/bi, and 
we then have 

b2fToK^ 


F ~ 


-ik 


(Cll) 


6i(s + 2Ao)™™ 

Leaving aside the constant factor in expr. (lC7p the corresponding anti-transformed fnnction 

is 




b2fT,K\.,_,. 


bi 


+0 - 


u 


exp [—2 Aom] 


(C12) 
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